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Abstract 



We present an analytic model for the interaction between planetary atmospheres and stel- 
lar winds from main sequence M stars, with the purpose of obtaining a quick test-model that 
estimates the timescale for total atmospheric mass loss due to this interaction. Planets in the hab- 
itable zone of M dwarfs may be tidally locked and may have weak magnetic fields, because of 
this we consider the extreme case of planets with no magnetic field. The model gives the plan- 
etary atmosphere mass loss rate as a function of the stellar wind and planetary properties (mass, 
atmospheric pressure and orbital distance) and an entrainment efficiency coefficient a. We use 
a mixing layer model to explore two different cases: a time-independent stellar mass loss and a 
stellar mass loss rate that decreases with time. For both cases we consider planetary masses within 
the range of 1 — > 10 M ffl and atmospheric pressures with values of 1, 5 and 10 atm. We apply 
our model to Venus by estimating its atmospheric mass loss rate by the interaction with the solar 
wind and compare our model with more detailed simulations. We find a good agreement between 
our results and the atmospheric mass loss obtained by more detailed models, and it is therefore 
appropriate for carrying out an exploration of the broad parameter space of exoplanetary systems. 
For the time dependent case, planets without magnetic field in the habitable zone of M dwarfs 
with initial stellar mass losses of < M w < 10~ n M Q yr _1 , may retain their atmospheres for at least 
1 Gyr. This case may be applied to early spectral type M dwarfs (earlier than M5). Studies have 
shown that late type M dwarfs (later than M5) may be active for long periods of time (> 4Gyr), 
and because of that our model with constant stellar mass loss rate may be more accurate. For these 
stars most planets may have lost their atmospheres in 1 Gyr or less because most of the late type 
M dwarfs are expected to be active. We emphasize that our model only considers planets without 
magnetic fields. Clearly we must expect a higher resistance to atmospheric erosion if we include 
the presence of a magnetic field. Nevertheless, as a first approximation our model is able to give 
a reliable time scale, as evidence by comparing our results with more detailed models. 

Keywords: Astrobiology; Atmospheres, evolution; Extrasolar planets; Solar wind 
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1 Introduction 



Low-mass main sequence stars of spectral class M (M dwarfs, or dM stars) are the most abundant stars 
in the galaxy, representing about 75% of the total stellar population. Their low masses (0.08 - 0.8 M©) 
allow them to have main sequence lifetimes of the order of 10 11 years (or more), i.e., significantly 
longer than the age o f the universe. Du ring this time the dM stars present nearly constant luminosities 
(L ~ 10 _1 - 10~ 3 L^. IScalo et al.ll2007|) . These peculiar characteristics make them interesting systems 
to search for habitable planets. However, there are some issues that make the habitability of planets 
around these stars a matter of debate. We will not re view here the habita bilit y issues rela t ed to p lanets 
around M stars, referring the reader to the papers by iTarter et al.1 (|2007|) and lScalo et al.1 (|2007|) . 



M stars exhibit significant temporal variability as a consequence of phenomena occurring in the 
region from their photospheres to their coronae. As a result, they emit large amounts of ultraviolet 
(UV) radiation and X-rays during their active periods and strong stellar winds. This stellar activity 
can have several effects on a planet. In this paper we focus on the interaction of stellar winds with 
the atmosphere of planets around these stars. Stellar winds are mainly charged particles (protons 
and electrons) ejected from the outer stellar atmosphere driven by the pressure-expansion of the hot 
corona. It has been calculated that the M dwarf stellar winds are denser and faster than winds from 



solar- like stars (IWargelin and Drake! . 120011) . and these winds may be able to erode the atmospheres of 
planets around the stars. 

The partial or total erosion of a planetary atmosphere affects directly the habitability of a planet, 
since a main requirement for the development of life on its surface is the presence of an atmosphere. 
The planetary atmosphere stabilizes the surface temperature and provides the pressure for water to be 
liquid, an essential requirement for life. Therefore, if a planetary atmosphere is removed by the stellar 
wind from its parent star (or any other process for that matter), this represents a serious problem for 
the habitability of the planet. For habitable planets around M dwarfs, the stellar wind may have a 
larger effect, not only because it is denser and faster than the winds from solar-type stars, but because 
their habitable zones (HZ) correspond to orbits very near to the star. 

The concept of HZ is commonly used in the context of planets that mi ght be suitable to k eep life 



on their surface. There are many definitions for the HZ (e.g lDoldl 19641 and lHeath et al.lll999|) but we 
take here Kasting's concept which is defined as the regi on around a s t ar in which a planet with an 
atmosphere is able to support liquid water on the surface (|Kasting et all 1 1993b . The formal definition 
is related to the stellar luminosity by: 



0.95AU < 



D 



< 131 AU 



(1) 



where D is the distance from the planet to the star, L t is the stellar luminosity and L Q is the solar 
luminosity. For the M dwarf case this distance is ~0.2 AU from the star. Planets in the HZ of M 



dwarfs would become tidally locked in about 1 Gyr or less (IGriefimeier et all 120001) . Calculations 
indicate that the magnetic moment of a tidally locked terrestrial planet may be 10 to 100 less intense 
than Earth's magnetic moment. Because it is likely that planets in the HZ of M dwarfs have weak 
magnetic fields, here we assume the extreme case where the planet has no intrinsic magnetic field, 
like Venus. 

The search for extrasolar planets is an interesting issue in the context of Astrobiology. The 



3 



smaller planets have been found orbiting dM stars, due to the fact that it is easier to observe ex- 
oplanets with the current techniques (gravitational lensing, radial velocities or transits) around less 
luminous stars. Since it is potentially possible to detect habitable planets around M dwarf stars, 
the study of the habitability conditions of these planets is decisive in the search for life around dM 
stars. Some recent exoplanet searches suc h as CoRoT (COnvection, Rotation & planetary Tran- 



sits) (lAgrain et all . [2008 



Barge et all 120081) and HA RPS (High Accuracy Radial velocity Planetary 



Search project) (lUdry et all 120071 : iPepe et all 12000) are dedicated to looking for planets with Earth- 



like masses around M dwarfs. 



2 Atmospheric mass loss calculation 

2.1 Erosion of the atmosphere by the stellar wind 

We first carry out an estimate of the atmospheric mass loss due to the interaction of a stellar wind of 
density p w and velocity v w with a planetary atmosphere. For an Earth-type planet, the atmosphere is 
a thin, gaseous layer covering the planet. Therefore, both the inner and outer radii of the atmosphere 
are « R P , where R P is the radius of the planet. 

For a planet with no magnetic field, the impinging stellar wind will form a bow shock around the 
body of the planet and its atmosphere. For low enough stellar wind densities this bow shock will 
not exist and the stellar wind will be directly absorbed by the planetary atmosphere. For a stellar 
wind/planet interaction in which the bow show does exist the post-bow shock material will have a 
sound speed c s ~ v w (as the shock will be strong due to the hypersonic nature of the stellar wind). The 
post-bow shock flow will envelop the leading hemisphere of the planet, and a turbulent mixing layer 
will be formed on the contact surface between the shocked wind and the planetary ionosphere. The 
shocked wind is stationary in the stagnation region of the bow shock, and has an increasing velocity 
towards the bow shock wings, reaching velocities of ~ v w towards the edge of the leading hemisphere 
of the planet/wind interaction flow. 



From the mixing layer formalism of ICanto and Ragal (|1991|) it is straightforward to show that a 



flow with a sound speed ~ v w and a tangential velocity of the same order has an entrained mass flux 
is M « ap w v w (where a is an entrainment efficiency which is determined through fitting the mixing 
layer model to laboratory experiments). Multiplying the entrained mass flux by the surface of the 
leading hemisphere, we then obtain an estimate for the atmospheric mass per unit time entrained by 
the stellar wind: 

M a ~ 2nR 2 P ap w v w (2) 

where M a represents the atmospheric mass loss. This equation can be combined with the relation 
M w = AnD 2 p w v w (where M w is the stellar mass loss rate and D is the orbital radius of the planet) to 
obtain: 



Through comparisons with laboratory experiments of plane, turbulent mixing layers. lCanto and Raga 



d 199 lh determined an a = 0.03 entrainment efficiency. It is of course not clear whether or not this 



value for a is appropriate for the different flow geometry of our stellar wind/planetary atmosphere 
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interaction model. Nevertheless. iBauer and Lammed (|2004) suggested that the value of a ~ 1/3 for 



the case of Venus. In order to sample the possible range of a values we use a =0.01 and 0.3. 

We should note that the interaction between the solar wind and Venus is found to be more complex 
than we have described above. For example, a magnetic barrier (of com pressed, magnetized solar 



wind material) is formed in the solar wind/planetary ionosphere interface (|Russell et all 120061) . This 
phenomenon could modify the properties of the entrainment process in a substantial way. Therefore, 
our entrainment model should be considered as a parametrization with which a range of possible 
entrainment rates can be explored regardless of the physical constraints for the a parameter. 

2.2 The timescale for total atmospheric erosion 

In order to estimate the initial atmospheric mass M atm as a function of the surface atmospheric pressure 
P , we use the hydrostatic equation for an isothermal atmosphere in the form: 

4nP Rl 

M atm = , (4) 

8 

where P is the surface pressure, R P the planetary radius and g the gravitational acceleration at the 
planetary surface. This equation can be combined with the relations g = GM P /Rp (where G is the 
gravitational constant and M P the planetary mass) and M P = 4nR p p P /3 (where p P is the average 
density of the planet) to obtain 

We now consider planets with different masses M P , assuming the Earth average density (i.e., p P = 
5.5g cm -3 ), and from Eq. © we obtain : 

M atm ( P Q \(M P xl 



Ma,* Uatm/UJ ' (6) 

where 1 atm=10 6 g cnT 1 s^ 2 , M ffi = 5.97 x 10 27 g and M a>fB = 5.27 x 10 21 g are the Earth's surface 
atmospheric pressure, planetary mass and atmospheric mass, respectively. 

From Eqs. ® and © we can then obtain the characteristic timescale t for total loss of the 
planetary atmosphere : 

M a llatm/\M P / V0.2AU/ \ a j\Mj * 

normalized to D = 0.2 AU, the average distance for a planet in the habitable zone around an M dwarf. 

If the mass loss rate of the star M w is time-independent, the timescale to is equal to the time for 
total erosion of the planetary atmosphere. From Eq. © we then conclude that the timescale for total 
atmospheric mass loss will satisfy the condition to > 1 Gyr provided that 

M w . I P n \(MA l/3 i D \ 2 /0.3 



M Q ~ 40 1 1 atm) \M P j V0.2AU/ [ a ' 
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Therefore, for an atmosphere with initial surface pressure P = 1 atm of a planet with M P = 1 M B and 
an orbital radius D = 0.2 AU, the timescale for total atmospheric loss will be long enough (« 1 Gyr) 
for the development of life if M w < 40 M Q , if we assume an a = 0.3 entrainment efficiency. Here M Q 
is the present solar mass loss, 10 -14 M Q yr _1 . For a lower, a = 0.01 entrainment efficiency (see §2.1), 
we derive a higher, M w < 1200 M Q value for the upper limit of the allowed stellar mass loss rate. 



2.3 Atmospheric loss due to a time-dependent stellar wind 



Observations indic ate that the stellar mass loss rate has a temporal dependence for solar-type stars 
(Wo od et all 120021) . We now assume a similar dependence on the loss rate for main sequence M stars. 



Wood et al.l (120021) obtained a temporal dependence following the power law 

2 



M w = M 



ir.O 



tw 



(9) 



wher e t w = 0. 1 Gyr and M W; o is the stellar mass loss rate at t = 0. For the Sun, M o = 2 x 10~ u M Q 
yr _1 ([Wood et all 120021) . In other words, in its early main sequence lifetime, the Sun had a mass loss 
rate ~ 2000 times larger than its present mass loss rate (10~ 14 Moyr^ 1 ). 
We now combine Eqs. © and © to obtain: 



Ma = 



dMa 
dt 



Rp 
D 



aM w , 



(t + t w ) 2 



(10) 



This equation can be integrated over time. Using Eq. © and the definition of Eq. © (t = M afi /M a ), 
we obtain: 



M a {t) = M afl , (11) 

\ to I t/t w + i 

where M a is the initial mass of the planetary atmosphere and t is given by Eq. © with M w = M w . 
By setting M a (t f ) = in Eq. (fTTT) we then obtain the time 



t/t M 



tolh 



1 - t /t w 



(12) 

t for t «: t w 



for total atmospheric mass loss. In Fig. 1 we plot tf as a function of t /t w . We have tf 
(= 0.1 Gyr, see above) and tf — » oo for t —* t w . 

From Fig. 1, we see that Eq. (TT21 gives tf/t w « 10 for to/t w « 0.9 (in other words, for to = 
0.9 x 0.1 Gyr). Therefore, tf > 1 Gyr (the timescale necessary for the development of life) is obtained 
for to > 0.09 Gyr. Using this limit in Eq. ©, we obtain the condition 



M 



vi'.O 



<433 



1 atm 



Me 
M P 



1/3 



D 



0.2 AU 



2 /0_3 
a 



(13) 



This condition allows for stellar wind mass loss rates a factor of ~ 10 higher than Eq. ®, which was 
derived assuming a time-independent stellar mass loss rate. 
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Figure 1: Total atmospheric erosion time tf resulting from the evolving stellar wind model (see Eq. 
fT2j) as a function of to (which is the total atmospheric loss time obtained assuming a time-independent 
stellar wind, see Eq. U}- Both tf and to are given in units of t w = 0.1 Gyr (see Eq. [9]). 
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3 Conditions for atmospheric survival obtained from the time- 
independent and time-dependent stellar wind models 



Figures 2, 3, 4, and 5 show the stellar mass loss necessary for a planet without magnetic field located 
at 0.2 AU to conserve its atmosphere for 1 Gyr, considering the time independent and time-dependent 
stellar wind models. Lower stellar mass loss rates than the ones given in these figures result in longer 
times for total loss of the planetary atmospheres (Eqs. [8]and[T3]). The mass loss rates shown in these 
figures have been computed from Eqs. © and (fT3b assuming the a = 0.01 entrainment coefficient (in 
the bottom of the range implied by laboratory experiments, see §2.1), the obtained mass loss rates are 
plotted in the Figs. 2 and 4. The results for a = 0.3 entrainment coefficient are presented in Figs. 3 
and 5. 

As expected, planets with larger initial atmospheric pressures have more resistance to erosion. 
More massive planets loose their atmospheres faster than smaller planets. This result is due to the fact 
that (for a set of planets with the same density p p , see Eq. 12.21) while the surface area presented to the 
stellar wind increases as M p 3 , the total atmospheric mass increases only as M p 3 (see Eq. [5]). 

From Figs. 2 to 5 we see that from the time-dependent stellar wind model (three upper lines) 
stellar mass loss rates required for totally loss a planetary atmosphere within 1 Gyr, are an order of 
magnitude higher than the ones obtained from the time-independent stellar wind model (three lower 
lines). These results imply that the atmosphere of a planet in the habitable zone of an active M dwarf 
may have an atmosphere that lasts for w lGyr only if the two following conditions are met : 

1. the entrainment parameter a is close to a = 0.01 

2. a substantial drop in the stellar wind mass loss rate occurs over a timescale of ~ 1 Gyr (in 
other words, that the mass loss rate falls at a rate comparable to the one assumed in our time- 
dependent stellar wind model, see §2.3). 

If these conditions are not met, only a planetary atmosphere with a very high initial atmospheric 
pressure may survive for timescales of ~ 1 Gyr under the action of a stellar wind with M w « 10~ 10 
M yr 1 . 



3.1 Venus comparison 



Several authors have suggested that Venus may have lost a considerable amount of its atmosphere 
trough its interaction with the solar wind. We used our time-dependent stellar wind model to calculate 
the atmospheric mass loss rate on the early Venus atmo sphere. The initial solar mass loss rate was 
considered to start with ~ 2 x 10~ 11 M Q yr~' at 0.1 Gyr (|Wood et all 120021) and evolve with time as 
presented in Eq. |9j We calculated the values for th e atmospheric mass l oss rate using Eq. \10\ 

Our results are comparable to those obtained by lKulikov et al.l(|2006[) . who obtained similar values 



for the atmospheric mass loss rate using a numerical test particle model (Fig. 6). iKulikov et al 
(|2006|) obtained series of values for the + pickup loss, considering the interaction of stellar winds 
and XUV flux of the young Sun, assuming different atmospheric density profiles. They used different 
parameters for the solar wind integrated over time for the period from 3.5 to 4.6 Gyr before present. 
They suggest that Venus may have lost during 4.5 Gyr more than 250 bars of + ions. In Fig. 6 we 
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Figure 2: Initial stellar mass loss rate needed in order to keep a planetary atmosphere for 1 Gyr in 
the habitable zone (at 0.2 UA) of a M dwarf, for planets with different masses. The three upper lines 
have been computed with the time-dependent wind model, and the three lower lines with the time- 
independent wind model. This figure corresponds to a = 0.01 and the solar mass loss considered here 
is for the current Sun value (M© « 10" 14 M© yr _1 ). 
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Figure 3: Initial stellar mass loss rate needed in order to keep a planetary atmosphere for 1 Gyr in 
the habitable zone (at 0.2 UA) of a M dwarf, for planets with different masses. The three upper lines 
have been computed with the time-dependent wind model, and the three lower lines with the time- 
independent wind model. This figure corresponds to a = 0.3 and the solar mass loss considered here 
is for the current Sun value (M© « 10~ 14 M Q yr~'). 
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Figure 4: Initial stellar mass loss rate needed in order to keep a planetary atmosphere for 1 Gyr in the 
habitable zone (at 0.2 UA) of a M dwarf, for planets with different atmospheric pressures. The three 
upper lines have been computed with the time-dependent stellar wind model, and the three lower lines 
with the time-independent wind model. This figure corresponds to a = 0.01 and the solar mass loss 
considered here is for the current Sun value (M© « 10" 14 M Q yr _1 ). 
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Figure 5: Stellar mass loss rate needed in order to keep a planetary atmosphere for 1 Gyr in the 
habitable zone (at 0.2 UA) of a M dwarf, for planets with different atmospheric pressures. The three 
upper lines have been computed with the time-dependent stellar wind model, and the three lower lines 
with the time-independent wind model. This figure corresponds to a = 0.3 and the solar mass loss 
considered here is for the current Sun value (M© « 10~ 14 M Q yr~'). 
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Figure 6: Comparative for the atmospheric mass loss rate on the early Venus atmosphere obtained by 
our time-dependent stellar wind model and more detailed models. The solid and dashed lines repre- 
sent the values corresponding to our mod el considering values of a=0.3 and 0.01 respectively. Squares 
and crosses show the values obtained by iKulikov et al.1 (|2006f) fo r their cases l b and 3 a, respectively 



Rhom bus, triangle and star rep resent the values calculated by iMoore et al.l (1199 1|) . iLammer et al 
(120061) and ILammer et al.1 (I2008|) . respectively. See text for more details. 



compare our model with cases lb and 3a from IKulikov et al.1 (120061) . Case lb is for the maximum 
solar wind and for the obstacle boundary that is one oxygen scale height below the exobase level, and 
case 3a is for the minimum solar wind and for the exobase level chosen as the boundary for the solar 
wind-exo sphere interaction. These cases were chosen because they have the maximum and minimum 
atmospheric mass loss rates for Venus. 

Our results were also compared with other works that estimate the pick up ra te for Q + ions consid- 
ering the present solar wind conditions. For example. ILammer et al.l (120060 and lLammer et al.l d2008|) 
use a gas dynamic test particle model to obta i n an a verage loss rate of 1.6xl0 25 s _1 and 1.5xl0 25 s~ l 
respectively. In an earlier work. IMoore et al.l (| 199 ll) estimated a value of about 1 x 10 24 s _1 applying 
a three-dimensional global hybrid simulation of the sola r wind interaction with the d ayside of Venus. 
Measure ments carried out by th e Pioneer Venus Orbiter dMihalov and Barneslll982|) and Venus Mars 



Express dBarabash et al.l . |2007|) indicate atmospheric mass losses of ~ 10 ions s ; therefore our 



model is in good agreement with more detailed models and observations of the Venus atmospheric 
mass loss, with the plot for a = 0.3 being closer to the results obtained by other authors. 
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3.2 Discussion: Survival of planetary atmospheres in the habitable zones of M 
dwarfs 



In gener al, there is little ob s ervati onal data for M dwarf winds. The first work on this topic was 
made by lWright and Barlowl (11975b . who pointed out that the stellar ma ss loss may be inferr ed from 
observations in the radio region of the spectrum. Following this idea, iMullan et al.l (| 19921) used a 
low frequency test with the James Clerk Maxwell Telescope (JCMT) and the Infrared Astronomical 
Satellite (IRAS). Their results suggested a stellar mass loss rate of a few times ~ 10 ~ 10 M ? yr" 1 , 
four orders of magnitude larger than the present solar mass loss. Later, ILim and White! (|1996l) used 
observations of the dMe stars YZ CMi and AD Leo at the millimetric wavelength range to show that 
the mass loss rate of this type of st ars must be less than ~ 10~ 13 M Q yr"' 



Van den Port and Dovldd 1997b analyzed the data of JCMT and IRAS reaching a different conclu- 



sion than lMullan et al.1 (|1992l) . By solving the radiative transfer equations for stellar winds from dMe 
stars, they showed that the inferred power-law flux distributions, based on radio, JCMT and IRAS 
data, cannot be reconciled with the flu x distributions from a stellar wind of - 
ing to 



10" 1U M Q yr . Accord- 

Van den Port and Dovle (119971) the maximum allowable mass loss rate is at most a few times 



~ 10" 12 M Q yr" 1 , while iHolzwarth and Jardind (120071) derive a maximum mass loss of 10M o , using a 
model for the wind properties of cool main -sequence stars, which com prises their wind ram pressures, 
mass fluxes, and terminal wind velocities J Wargelin and Drake d200ll) proposed a new observational 
technique for measuring stellar winds from late type dwarfs (F-M) through o bservations of charge 
exchange-induced X-ray em ission, but its results have been limited until now (|Wargelin and Drakd . 



2002l : IWargelin et all 120081) . 



Many M dwarfs present their Ha line in emission and because of that are classified as dMe stars. 
The strength of the Ha line is considered as a proxy of activity associated to higher mass loss rates 
than those for stars with no Ha emission. Pne of th e most observed dMe stars i s Proxima Centauri that 
has a measured mass loss rate < 3 x 10" 13 M Q yr" 1 (|Wargelin and Drakd . 12002). Proxima Centauri is a 
M5.5Ve star member of a multiple system with an age of 4 - 4.5 x 10 9 years (IDemarque et all 1 19860 . 
Using this range for the age of Proxima Centauri and Eq. ©, the initial mass loss rate of Proxima 
Centauri lies between 5x l0~ 1Q M n yr" 1 and 6.4 xl0~ 10 M Q yr" 1 . If Proxima Centauri is a young star 
(< 1 Gyr) as obtained by Fleming et al.1 (| 1995b . then its initial mass loss rate would be ~ 4 x 10 -11 
M Q yr" 1 . For the case of other dMe sta rs that are considered young, like YZ CMi (dM4.5e) and 



AD Leo (dM3.5e) (Fleming et all 11995), the calculated maximum for their mass loss rate is 10 



-12 



M yr (Van den Port and Doy Id 11997b. The ag e of these stars lies somewhere between 0.5 Gyr 



dSoderblom L 1 1990b and 1 Gyr (lYoung et all 119890 . Therefore their initial mass loss rates may be in 
the range of l.lxlO" 11 M Q yr" 1 to 3.6xl0" n M Q yr" 1 . YZ CMi and AD Leo are among the most 
active M dwarfs, and we would therefore expect that they represent an upper limit for M dwarf mass 
loss rates. 

Here we are considering that the mass loss rate for dM stars is time depende nt, but obse r vation s 
indicate that some M dwarfs may not follow the time dependence proposed by I Wood et al.l (|2002h . 
For example. ISilvestri et al.l (|2007b found that the lower mass M dwarfs (spectral types M3 and later) 
are more likely to be active at old ages than higher mass M dwarfs (M0-M3V M stars in clusters 



with spectral type s M3 and later seem to remain active for more than 4 Gyr (ISilvestri et al.l . |2007J) 



Gizis et al. (2000) found that late type M dwarfs (later than M6.5) that are strong Ha emitters are 
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likely to be old (>1 Gyr). Then, for the analysis of our results we will consider that the time dependent 
mass loss rate applies to earlier M dwarfs (earlier than M5) and the constant mass loss rate fits better 
the behavior of the later spectral type M dwarfs (later than M5). For t hese late type stars the percentage 
of act ive stars goes from ~40% for M5 stars to ~60% for M7 stars (|Silvestri et all 120071 : iGizis et al. . 



2000). 



For those early type M dwarfs, the atmosphere of a planet in their HZ may be safe for stellar mass 
loss rates < 5 x 10 3 M o = 5 x 10~ n M Q yr" 1 for a = 0.01 (Figs. 2 and 4). For a = 0.3, only planetary 
atmospheres with more than 5 atm may be able to survive for 1 Gyr or more if the parent star of the 
planet has mass losses < 1 x 10 3 M o = 1 x 10~ n M yr" 1 (Figs. 3 and 5). Because most of early type 
M dwarfs are not active (|Silvestri et all 120071 : IGizis et alll2000r) we would expect that their mass loss 
rates would be less than the observed ~ 10 -11 M Q yr 1 for the most active M dwarfs (YZ CMi and AD 
Leo). We therefore expect that most planets in the habitable zone of early type M dwarfs may keep 
their atmospheres for 1 Gyr or more. 

For M dwarfs that do not have a mass loss rate that decreases as a function of time our time- 
independent wind model predicts the worst scenario for planets in the habitable zone of these stars. 
If a = 0.01 a planet must have an atmosphere of 5 atm or more in order in order to survive mass loss 
rates of 10~ n M Q yr^ 1 (Figs. 2 and 4). For a = 0.3 (regardless of the planetary mass/atmospheric 
mass values) a planet in the HZ of an active late type M dwarf will loose its atmosphere in 1 Gyr or 
less (Figs. 3 and 5). 



4 Conclusions 

Because the habitable zone of M dwarfs is very close to the star (<2 AU) planets in that HZ may be 
tidally locked and therefore sl ow rotators. As a consequence these planets may have weak magnetic 



fields (IGrieBmeier et all 120001) . In this work we have considered the extreme case of planets with no 
magnetic field. 

Our model gives the planetary atmosphere mass loss rate as a function of the stellar wind and 
planetary properties (mass, atmospheric pressure and orbital distance) and an entrainment efficiency 
coefficient a. The model is applied to explore two different cases: a time-independent stellar mass 
loss and a stellar mass loss rate that decreases with time. We find that our time-independet stellar 
wind model predicts similar results to the ones found by using detailed models of the atmospheric 
erosion of Venus. Therefore, our model is valid to test if a planet is able to keep its atmosphere for 
a certain period of time. In other words, if one planet (under the assumption that the planet has no 
magnetic field) doesn't pass our test, it is not necessary to analyze the atmospheric erosion with more 
detailed models. 

Based on observations of M dwarfs we derive that our time dependent model may better describe 
early type M dwarfs (earlier than M5) and the case for constant stellar mass loss more appropriate 
for the interaction of M dwarfs with spectral types later than M5 with the atmospheres of planets 
in their habitable zones. Our results indicate that early type M dwarfs may are more promising for 
habitable planets searches because there are several planetary mass/atmospheric mass combinations 
for which a planet without magnetic field in the habitable zone may preserve its atmosphere for at 
least 1 Gyr. Planets in the habitable zone of active late type stars (later than M5) are less likely to 
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keep their atmospheres for 1 Gyr or more. 

It is clear that more observations are needed to narrow down the mass loss rate values for M 
dwarfs. Particularly, such mass loss rate determinations will have to be made for M dwarfs in which 
the presence of planets is indeed detected, so as to be able to determine (from our model) whether or 
not the planets are candidates for the presence of life. 

We end the discussion by noting that our model for the erosion of planetary atmospheres is re- 
stricted to the case of planets with no magnetic fields. Clearly, the presence of a planetary magneto- 
sphere would protect the atmosphere, and may allow the planet to preserve its atmosphere for longer 
timescales in the presence of stellar winds with higher mass loss rates. The study of the survival of 
the atmosphere of a magnetized planet under the interaction with an M dwarf wind is left for a future 
study. 
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